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Abstract 



. We use gravitational lagrangians ROi^ R^/ —g and linear combina- 

tions of them; we ask under which circumstances the de Sitter space- 
time represents an attractor solution in the set of spatially flat Fried- 
man models. 

Results are: For arbitrary k, i.e., for arbitrarily large order 2k + 4 
of the field equation, on can always find examples where the attractor 
property takes place. Such examples necessarily need a non-vanishing 
i?^-term. The main formulas do not depend on the dimension, so one 
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gets similar results also for 1+1-dimensional gravity and for Kaluza- 
Klein cosmology. 

PACS number: 0450 Unified field theories and other theories of gravita- 
tion 

1 Introduction 

Over the years, the notion "no hair conjecture" drifted to "no hair theo- 
rem" without possessing a generally accepted formulation or even a complete 
proof. Several trials have been made to formulate and prove it at least for 
certain special cases. They all have the overall structure: "For a geometri- 
cally defined class of space-times and physically motivated properties of the 
energy-momentum tensor, all the solutions of the gravitational field equation 
tend asymptotically to a space of constant curvature." 

1.1 Historical notes 

It is the aim of this paper to clarify the relations between the several exist- 
ing versions, and then to develop the cosmological no hair theorem towards 
applicability to a certain class of higher order field equations. Let us start 
with some historical notes. 

The paper [1] by Weyl (1927) is cited in [2] with the phrase "The be- 
haviour of every world satisfying certain natural homogeneity conditions in 
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the large follows the de Sitter solution asymptotically." to be the first pub- 
lished version of the no hair conjecture. 

Barrow and Gotz [2] apply the formulation "All ever-expanding universes 
with A > approach the de Sitter space-time locally." 

(Ever-expanding to be meant as: there is a time such that for all 
t > to the Hubble parameter is positive. In other words: a bounce is allowed, 
a recoUapse is not allowed.) 

Let us comment this formulation: So they circumvent the necessity to 
distinguish the initial data between expanding and recoUapsing ones, but 
their formulation needs a further explanation; example: If one changes the 
initial data continuously from recollapsing to ever-expanding ones, then one 
gets a critical value of the initial data between them, where one has an ever- 
expanding model which need not tend to the de Sitter space-time but can 
have typically a linear expansion law. So these critical values of the initial 
data have to be excluded, too. 

The first proof of the stability of the de Sitter solution (here: within 
the steady-state theory) , is due to Hoyle and Narhkar [3] . In the papers [4] 
by Price perturbations of scalar fields have been considered for the no hair 
theorem. The probability of inflation is large if the no hair theorem is valid, 
cf. [5]. 
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1.2 Physical properties 

Peter, Polarski and Starobinsky [6] compared the double-inflationary mod- 
els with cosmological observations. Barrow and Gotz discussed the no hair 
conjecture within Newtonian cosmological models [2,7]. Hiibner and Ehlers 
considered inflation in an open Friedman universe and have noted that in- 
flationary models need not to be spatially flat [8]. 

Gibbons and Hawking have found two of the earliest strict results on 
the no hair conjecture for Einstein's theory [9] in 1977. Barrow gave ex- 
amples that the no hair conjecture fails if the energy condition is relaxed 
and points out, that this is necessary to solve the graceful exit problem. He 
uses the formulation of the no hair conjecture "in the presence of an effec- 
tive cosmological constant (e.g. from viscosity) the de Sitter space-time is a 
stable asymptotic solution" . This is a much weaker statement because only 
space-times in a neighbourhood of the de Sitter space-time are involved. He 
mentioned that an ideal fluid with equation of state p = —p is equivalent to 
a A-term in some cases but not always [10]. 

Usually, energy inequalities are presumed for formulating the no hair 
conjecture. Nakao, Shiromizu and Maeda [11] found some cases where it 
remains valid also for negative Abbott-Deser mass [the latter goes over to 
the well-known ADM- mass (Arnowitt, Deser, Misner) for A — > 0]. They 
cite Murphy [12]. In [12], viscosity terms as source are considered to get a 
singularity-free cosmological model. Murphy [12] used Einstein's theory, and 
Oleak [12] made similar considerations within Treder's theory of gravity. 
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In the Eighties, these non-singular models with viscosity where re-inter- 
preted as inflationary ones, cf. e.g. [13]. 

In the three papers [14] Prigogine et al. developed a phenomenological 
model of particle and entropy creation. It allows particle creation from space- 
time curvature, but the inverse procedure (i.e. particle decay into space-time 
curvature) is forbidden. This breaks the t — > — i-invariance of the model. 
Within that model, the expanding de Sitter space-time is an attractor solu- 
tion independently of the initial fluctuations; this means, only the expanding 
de Sitter solution is thermodynamically possible. To these papers cf. also 
[15]. 

Vilenkin [16] discussed future- eternal mRating universe models; they must 
have a singularity if the condition D: "There is at least one point p such that 
for some point q to the future of p the volume of the difference of the pasts 
of p and q is finite" is fulfilled. 

Mondaini and Vilar have considered recollapse and the no hair conjecture 
in closed higher-dimensional Friedman models [17]. PuUin [18] discussed re- 
lations between the onset of black hole formation and the no hair conjecture. 
Concerning the no hair conjecture Shiromizu, Nakao, Kodama and Maeda 
[19] gave the following argument: If the matter distribution is too clumpy, 
then a large number of small black holes appears. Then one should look for 
an inflationary scenario where these black holes are harmless. They cannot 
clump together to one giant black hole because of the exponential expansion 
of the universe; this explains the existing upper bound of black holes in the 
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quasi-de Sitter model. Shibata, Nakao, Nakamura and Maeda have consid- 
ered asymptotic gravitational waves in an axially symmetric quasi de Sitter 
space-time [20] . They use numerical methods. The magnitude of black holes 
is restricted: above Merit — there do not exist horizons; this restriction 
one gets by considering a perturbed Schwarzschild-de Sitter-solution. The 
cosmic hoop conjecture expresses that the mass of a black hole in a quasi 
de Sitter model is bounded from above by Merit — and its surface is 
analogously restricted. 

The notion "quantum hair" means quantum numbers presenting quantum 
fields which should be classically forbidden if the no hair theorem is vahd. 
Coleman, Preskill and Wilczek found examples of quantum hairs on black 
holes [21]. 

Xu, Li and Liu [22] proved the instability of the anti-de Sitter space-time 
(classical instability against gravitational waves, and dust matter perturba- 
tions); one has A < 0, and in an open Friedman model the scale factor a 
in dependence of synchronized time t reads a = a cos ^ where A = — 3/q;^. 
The anti-de Sitter model has closed time-hke curves everywhere; a Cauchy 
horizon is the surface where closed time-like curves begin to exist, and there- 
fore, the anti-de Sitter model has no Cauchy horizon. (Of course, a closed 
curve has no beginning; the formulation means: The Cauchy horizon is the 
topological boundary of the set of point possessing the property that they 
are contained in a closed time-like curve.) 

Coley and Tavakol discussed the robustness of the cosmic no hair conjec- 
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ture under using the concept of the structural stabihty [23] (compare with 
chapter 5 below). 

1.3 Fourth-order gravity 

Sirousse-Zia considered the Bianchi type IX model in Einstein's theory with a 
positive A-term and got an asymptotic isotropization of the mixmaster model 
[24]. She cites (and uses methods of) Belinsky, Lifshitz and Khalatnikov [25]. 
Miiller [26] used L = and discussed the power-asymptotes of Bianchi 
models. Barrow and Sirousse-Zia [27] discussed the mixmaster i?^-model 
and the question, under which conditions the Bianchi type IX 
model becomes asymptotic de Sitter ? 

Yokoyama and Maeda [28] considered the no hair conjecture for Bianchi 
type IX models and Einstein's theory with a cosmological term. They dis- 
cussed B? inflation in anisotropic universe models and got as a 

result that typically, an initial anisotropy helps to enhance inflation. For 
Bianchi type IX they got some recoUapsing solutions besides those converging 
to the de Sitter solution. 

Cotsakis, Demaret and de Rop [29] discussed the mixmaster universe in 
fourth-order gravity. To take the metric diagonal they write "is probably 
well justified" ; they discuss all types of curvature-squared terms. Paper [27] 
is continued in [30] by Spindel, where also general Bianchi type I models in 
general dimensions are considered. 

Gurovich et al. [31] considered L — R -\- aR^/^ to get a singularity- 
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free model in 1970, the solutions are of a quasi de Sitter type. One should 
remember that in spite of de Sitter's papers in the twenties, the inflationary 
cosmological model became generally accepted in only in 1979/80. 

The papers [32] , [33] consider the no hair conjecture for E? models, they 
use the formulation " asymptotical de Sitter, at least on patch" . The restric- 
tion "on patch" is not strictly defined but refers to a kind of local vahdity 
of the statement, e.g., in a region being covered by one single synchronized 
system of reference in which the spatial curvature is non-positive and the en- 
ergy conditions are fulfilled. The Starobinsky model is outlined as one which 
does not need an additional inflaton field to get the desired quasi de Sitter 
stage. One should observe a notational change: There, L = R + aR^ In R 
was called Starobinsky model, whereas L — R-\-aR? got the name "improved 
Starobinsky model" - but now the latter carries simply the name " Starobin- 
sky model". (For the inflationary phase, both versions are quite similar.) 
A further result of the papers [32] is that by the addition of a cosmologi- 
cal term, the Starobinsky model leads naturally to double inflation. Let us 
comment this result: It is correct, but one should add that this is got at 
the price of getting a "graceful exit problem" (by this phrase there is ment 
the problem of how to finish the inflationary phase dynamically) - in the 
Starobinsky model this problem is automatically solved by the fact that the 
quasi de Sitter phase is a transient attractor only. The papers [34] discuss 
the no hair conjecture within i?^-models and found inflation as a transient 
attractor in fourth order gravity. The papers [32] and [35], [37] discuss the 
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stability of inflation in i?^-gravity. The papers [36] discuss generalized cosmic 
no hair theorems for quasi exponential expansion. In Starobinsky [38] the 
no hair theorem for Einstein's theory with a positive A-term is tackled by 
using a sequence as ansatz to describe a general space-time. However, the 
convergence of the sequence is not rigorously proven. 

Starobinsky and Schmidt [39] have generahzed the ansatz of Starobinsky 
[38] to consider also the no hair theorem for L = R^. 

Shiromizu et al. [19] discussed an inflationary inhomogeneous scenario 
and mention the open problem how to define asymptotical de Sitter space- 
times. In Pacher [40] it is mentioned that only a local version of this conjec- 
ture can be expected to hold true, and that neither the definition of asymp- 
totic de Sitter nor the necessary presumptions to the energy-momentum ten- 
sor are clarified - two problems which are not finally solved up to now. The 
authors of [41] consider the no hair theorem for a special class of inhomoge- 
neous models and give partial proofs. Morris [42] considers inhomogeneous 
models for R + i?^- cosmology. In [43] inflation in inhomogeneous but spheri- 
cally symmetric cosmological models is obtained only if the Cauchy data are 
homogeneous over several horizon lengths. The analogous problem is consid- 
ered in [44] also with inclusion of colliding plane gravitational waves, they 
give a numerical support of the no hair conjecture by concentrating on the 
dynamics of gravitational waves. 

Berkin [45] gets as further result, that for L = f{R), a diagonal Bianchi 
metric is always possible. Similarly, Barrow and Sirousse-Zia [27] and Spin- 
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del [30] worked on diagonalization problem. They apply the diagonalisation 
condition of MacCallum et al. [46]. In 1918 Kottler [47] found a simple 
closed-form static spherically symmetric vacuum solution for Einstein's the- 
ory with A-term in Schwarzschild coordinates. 

ds^ = A{r)dt^ - - r\de^ + sin" 9 d^^) (1.1) 

with A{r) = 1 — "^ — ^r^. At the horizon /I = the Killing vector changes 
its sign and one gets by interchanging the coordinates t and r the correspond- 
ing Kantowski- Sachs model. Moniz [48] (1993) discusses the cosmic no hair 
conjecture within Kantowski-Sachs models and A > 0. He gets the de Sitter 
space-time not only asymptotically, but exactly in an anisotropic 3-|-l-slicing 
of space-time. He discusses the initial data that lead to a recoUapse and 
find them to be very rare; but the measure he uses is not well-defined, so, 
possibly, this is not the last word. It is curious to observe that he works 
with complicated elliptic integrals instead of applying the Schwarzschild-de 
Sitter-solution found by Kottler [47] in 1918. 

[49] gives an overview about the geometry of the de Sitter space-time. 

1.4 Sixth and higher order models 

The paper [50] by Buchdahl (1951) deals with lagrangians of arbitrarily high 
order. Its results are apphed in paper [51] to general Lagrangians F{R, □). 
From another motivation, BoUini et al. [52] consider higher-order field theo- 



10 



ries of the type 

and give solutions in the sense of distributions. 

Forgacs et al. [53] consider the non-local lagrangian R^R — M as Wess 
Zumino Witten model. 

The paper [54] by Vilkovisky was presented at the A. Sacharov-memorial 
conference held in Moscow in May 1991. In [54], the Sacharov-approach was 
generalized. The original idea of Sacharov (in 1967) was to define higher 
order curvature corrections to the Einstein action to get a kind of elasticity 
of the vacuum. Then the usual breakdown of measurements at the Planck 
length (such a short de Broglie wave length corresponds to such a large mass 
which makes the measuring apparatus to a black hole) is softened. Vilkovisky 
discusses the effective gravitational action in the form /?/(□)/?, where 

/(□) ^ J ^p{x)dx 

Martin and MazziteUi [55] discuss the non-local Lagrangian R^R as con- 
formal anomaly in two dimensions. 

Let us now come to sixth-order equations. Stelle [56] (1977) considers 
mainly fourth order i?^-models; in the introduction he mentioned that in the 
next order, terms like R^ + Rij-^R^^''' become admissible, but the pure R^- 
term is not admissible. Treder [57] used higher-order lagrangians, erspecially 
i?^-terms, and he mentioned that for R+R^iR ^g^^ a sixth-order field equation 
appears. Remark: This lagrangian leads to the same field equation as 
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R - ROR. 

In [58] , inflationary models with a term R^iR'^ in the action are considered, 
but they do not vary with respect to the metric, and so no sixth-order term 
in the field equation appears. 

Lu and Wise [59] consider the gravitational Lagrangian as a sequence 
S — S0 + S1 + S2 + . . . ordered with respect to physical dimension. So, 5*0 = -R 
and 5*1 sums up the R^-terms. They try to classify the 5'2-terms; however, 
their identity (8) is not correct, so they erroneously cancel the essential term 
RnR. 

Kirsten et al. [60] consider the effective lagrangian for self-interacting 
scalar fields; in the renormalized action, the term 

_DR_ 
c + R 

appears. Wands [61] classifies lagrangians of the type F{R, and men- 

tions that not all of them can be conformally transformed to Einstein's the- 
ory. Ref. [62] considers the lagrangian $^ni?, [63] the Lagrangian ROR^ [64, 
65] double inflation from $ and i?^ -terms, also the ROR-terms is discussed. 
Besides ROR Berkin [65] considers the de Sitter space-time as attractor solu- 
tion for field equations where the variational derivative of the term CijkiC'^-'^^ 
is included. The state of the art of the lagrangian ROR can be found in the 
papers [64-68]. 

The paper is organized as follows: Set. 2 compares several possible defini- 
tions of an asymptotic de Sitter space-time, set. 2.1. for the set of spatially 
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flat Friedman models, set. 2.2. for less symmetric models. Set. 3 deals with 
the Lagrangian and corresponding field equations for higher-order gravity. 
In set. 4, we determine under which circumstances the Bianchi models in 
higher-order gravity can be written in diagonal form without loss of gen- 
erality; the answer will be more involved than the analogous problem for 
General Relativity. In set. 5 we discuss the results from the point of view of 
structural stability in the sense of the "Fragility" -paper [23]. 

2 Definitions of an asymptotic de Sitter 
space— time 

In this section we want to compare some possible definitions of an asymptotic 
de Sitter space-time. 

2.1 Spatially flat Friedman models 

Let us consider the metric 



which can be called spatially flat Friedman model in n spatial dimensions. 
We consider all values n > 1, but then 

concentrate on the usual case n = 3. If n < 3 we often write a;, y, and 
z instead of x^, x^, and x^, resp. For metric (2.1) we deflne the Hubble 



n 





i=l 
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parameter H — a = ^. We get 



Roo^-n{— + H^), R^-2n{—+mH^) 
at at 



where m :— > 1. We get 
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Lemma 1: The following conditions for metric (2.1) are equivalent. 
A: It is flat. B: R = Rqo = 0. C: RijR'^ = 0. 
H: a — const, or [n = 1 and a — ln\t — to\ + const. ] 

Proof: A =^ B is trivial; B ^ D is done by solving the corresponding 
differential equation; D =^ A is trivial for a — const., the other case, i.e., 
ds'^ — dt^ — {t — toYdx^, represents flat space-time in polar coordinates; 
C B follows from the identity 

R,jR'^ = (i?oo)' + -{R - Roof (2.3) 

An analogous statement can be formulated for the de Sitter space-time. 
It holds 

Lemma 2: The following conditions for metric (2.1) are equivalent. 
A: It is a non-flat space-time of constant curvature. 

B: i?oo = 4t = ^'^^^^^ ^ 0- 

C: (n + l)RijR'^ = i?^ = canst. ^ 0. 

D: = const. 7^ or [n = 1 and ds^ = dt^ — &\v?{Xt)dx^ or ds^ = dt^ — 
smh^{Xt)dx^ ] 

The proof is analogous to lemma 1. 

For n — 1, the de Sitter space-time and anti-de Sitter space-time differ by 
the factor (-1) in front of the metric only. For n > 1, under the presumption of 
lemma 2, only the de Sitter space-time {R < 0) is covered. Lemma 2 shows 
that within the class of spatially flat Friedman models, a characterization 
of the de Sitter space-time using polynomial curvature invariants only, is 
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possible. 

Next, let us look for isometries leaving the form of the metric (2.1) in- 
variant. The function 

a{t) — c + a{±t + to), c, to — constants, (2-4) 

leads to an isometric space-time. The simplest expressions being invariant by 
such a transformation are and H. We take a as dimensionless, then H is 
an inverse time and H an inverse time squared. Let if 7^ in the following. 
The expression 

s HH-^ (2.5) 

is the simplest dimensionless quantity defined for the spatially flat Friedman 
models (2.1) and being invariant with respect to the isometries (2.4). Let 
n > 1 in the following: Two metrics of type (2.1) are isometric if and only 
if the corresponding functions a and a are related by equation (2.4). All 
dimensionless invariants containing at most second order derivatives of the 
metric can be expressed as f{s), where / is any given function. But if one 
has no restriction to the order, one gets a sequence of further invariants 

e, = HH-\...,e, = ^H-P-' (2.6) 

It holds: Metric (2.1) with H ^ represents the de Sitter space-time iff 
£ = 0. 

A third possible approach is the following: a{t) = Ht 

with H — const. 7^ is the de Sitter space-time, so we define an asymp- 
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totic de Sitter space-time by the condition 

lim ^ = const. ^ (2.7) 

Let us summarize the variants Var(i) of the definitions. 

Definition: Let i7 > in metric (2.1) with n > 1. We call it an asymptotic 

de Sitter space-time if 

Var (1): lim^^oo = const. > 

Var (2): limt-^oo — const. > and limt_^oo(?T' + l)RijRij — R^ — 
Var (3p): for 1 < j < p it holds hm^^oo^p = 0. All these definitions are 
different. One uses RijR^^ — n^{H -\- H'^Y -\- n{H + nH'^y. However, as 
we will see in section 3.2, all these definitions lead to the same result if we 
restrict ourselves to the set of solutions of the higher-order field equations. 

2.2 Inhomogeneous cosmological models 

Let us start looking at the Kottler metric eq. (1.1). The critical mass Merit = 
mentioned in subsection 1.3. in connection with the hoop conjecture can 
be deduced (at least for the symmetries of the Kottler metric as follows: At 
a horizon, the function A 

must vanish. One can see from eq. (1.1) that for A > the equation 
^4 = has solutions with positive values r if and only if m < Merit- This 
means: the hoop conjecture is valid in the class of spherically symmetric 
solutions. 

However, this is not the problem we are dealing with here. The 
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problem is that none of the above definitions can be generahzed to inho- 
mogeneous models. One should find a polynomial curvature invariant which 
equals a positive constant if and only if the space-time is locally the de Sitter 
space-time. To our knowledge, such an invariant cannot be found in the lit- 
erature, but also the non-existence of such an invariant has not been proven 
up to now. 

This situation is quite different for the positive definite case: For signature 
(+ + ++) and Sij = Rij - ^gij it holds: 

I={R- Rof + CijkiC'^"' + SijS'^ = 

iff the V4 is a space of constant curvature i?o- So / — > is a suitable 
definition of an asymptotic space of constant curvature. 

One possibility exists, however, for the Lorentz signature case, if one al- 
lows additional structure as follows: An ideal fiuid has an energy-momentum 
tensor 

Tij = {p + p)UiUj - pQij 

where Ui is a continuous vector field with UiU^ = 1. For stiff matter (p = — p), 
the equation T*-*.^- = implies p — const., and so every solution of Einstein's 
theory with stiff matter is isometric to a vacuum solution of Einstein's theory 
with a cosmological term. The inverse statement, however, is valid only 
locally: 

Given a vacuum solution of Einstein's theory with a A-term, one has to 
find continuous time-like unit vector fields which need not to exist from topo- 
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logical reasons. And if they exist, they are not at all unique. So, it becomes 
possible to define an invariant J which vanishes iff the space-time is de Sitter 
by transvecting the curvature tensor with u^u^ and/or g'^^ and suitable linear 
and quadratic combinations of such terms. Then time t becomes defined 
by the streamlines of the vector u^. If one defines the asymptotic de Sitter 
space-time by J — > as t — > oo, then it turns out, that this definition is 
not independent of the vector field u^. 

3 Lagrangian F{R, UR^ U^R^ . . . , U^R) 

Let us consider the Lagrangian density L given by 

L = F(i?, Di?, □2i?,...,n^'i?)y=^ (3.1) 

where R is the curvature scalar, □ the D'Alembertian and gij the metric 
of a (Pseudo-)Riemannian Vd of dimension D > 2 and arbitrary signature; 
g — — \ det gij \ . The main application will he D = A and metric signature 
(H ). F is supposed to be a sufficiently smooth function of its argu- 
ments, preferably a polynomial. Buchdahl [50] already dealt with such kind 
of Lagrangians 

in 1951, but then it became quiet of them for decades. Since 1990 a 
sequence of papers on this topic appeared: refs. [51, 68] for general k, and 
refs. [53 - 67] for the special case k — 1, i.e. the Lagrangian scalar is 
F{R, OR). 
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3.1 The field equation 

The variational derivative of L with respect to the metric yields the tensor 

P'' = 1^ (3-2) 

The components of this tensor were given in the first paper of ref. [51], their 
covariant components read 

1 

Pij = GRij — ^FQij ~ G-ij + gijOG + Xij (3.3) 

where the semi-colon denotes the covariant derivative, Rij the Ricci tensor, 
and 

X^J = E l9^AFA{B^-'Rrn■,m ' Fj,^,,P^-' R],,) (3-4) 

having the round symmetrization brackets in its last term. For /c = 0, i.e. 
F — F{R), a case considered in set. 4, the tensor X^j identically vanishes. It 
remains to define the expressions Fa, A — 0, . . . ,k . The definition given in 
[51] can be simplified as follows 

OF 

Fj^ = -TL— (3.5) 

and for A = A; — 1, . . . , 

dF 

Fa = DFa+i + 7—^ (3.6) 
and finally G — Fq. The brackets are essential, for any scalar $ it holds 

□($,) - (□$), = (3.7) 
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Inserting $ = D'^R into this equation, one gets identities to be applied in 
the sequel without further notice. 
It is well-known that 

P);i = (3.8) 

and Pij identically vanishes if and only if F is a divergence, i.e., locally 
there can be found a vector such that F = v^.^ holds. (Remark: Even for 
compact manifolds without boundary the restriction " locally" is unavoidable; 
example: Let D = 2 and V2 be the Riemannian two-sphere 5*^ with arbitrary 
positive definite metric, i? is a divergence, but there do not exist continuous 
vector fields fulfilling R = on the whole S'^.) 
Example: for m, n > it holds 

O^'R □"i? - R = divergence. (3.9) 

So, the terms □"^i? n'^R with naturals m and n can be restricted to the case 
m = without loss of generality. However, the more far-reaching statement 
by Wands [61, page 271] "Thus I can take any polynomial F{n^R) to be lin- 
ear in its highest-order derivative □"i?, multiphed by Fn{R)" is not correct. 
Let us give a counterexample: RDRnR, which leads to an eighth-order field 
equation. 

Proof that this is a counterexample: From dimensional reasons only ingredi- 
ents with < length >~^° are to be considered. Neglecting the divergencies, 
only the following ones are candidates: R^, R^OR, R^U^R, RU^R. They 
give rise to field equations of orders 4, 6, 8, and 10 resp. So the last term 
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cannot be included. It remains to look for 

F = RORDR + -fR^ + PR^DR + aR^U^R 

The variation of F with respect to the metric should vanish identically. Van- 
ishing of the 8th-order term requires a = — |. Vanishing of the 6th-order 
terms gives rise to the equation 

(□i? + ^i?^)(nit:),, = o 

For no value of (3 this is identically satisfied. 
3.2 Higher-order gravity 

We will examine the attractor property of the de Sitter space-time in the set 
of the spatially flat Friedman models. We need some useful relations for the 
de Sitter space-time: 

i? = + 1)//^ (3.10) 

and 

R 

Rij = —^9ij and d'^R = for A; > . (3.11) 
We insert this into the the field equation (3.3) 

= GR,j - ^Fg,j = g,, (^^^^ " (3-12) 

for the de Sitter space-time. The de Sitter space-time solves the field equa- 
tion if and only if 2RG ~ DF. If we choose the lagrangian {—R)^ with 
w e R the £)-dimensional de Sitter space-time satisfies the field equation 
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iS u— = ^. We will examine the attractor property of the de Sitter 
space-time in the set of the Friedman models for the lagrangian (— i?)" with 
2u — D — n + l>2. From this lagrangian it follows 



(3.13) 



and 



We get the field equation 



G = -u{-R) . 



(3.14) 



= -u{-Rr^'R,,--g,j{-Rr+u 



-R) 



M-l 



g,,un {-RY^' . (3.15) 



It is enough to examine the 00-component of the field equation, because all 
the other components are fulfilled, if the 00-component is fulfilled. We make 

the ansatz 

a{t) = 1 + P{t) (3.16) 

and get 



Rqq — —n$ — 2n(3 — n 
R = -2nP-2{n^ + n)(3-{n^ + n) (3.17) 
{-R)"^ = 2nm{n'^ + n)"^-^/3 + 2m{n'^ + n)"*/3 + (n^ + n)"* . 

One gets the field equation 

= -2n'^u{u + l)(n^ + n)""^/3 - 2n\{u - l)(n^ + n)u -2(3 + 

+nM(2M-n-l)(n2+n)"-V+(nK-^(n2+n)^(n2+n)"-^ (3.18) 
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Using the condition 2u — n + 1 we get 

= /3 + n/3. (3.19) 
All solutions of the linearized field equation are 

p{t) = Ci + 026-""' . (3.20) 

It follows 

a{t) ^t + cit + C2e-"* + C3 (3.21) 

and 

lim ^ = 1 + ci . (3.22) 

t—^00 t 

The D-dimensional de Sitter space-time is an attractor solution for the la- 
grangian F = (—R)^. 

The lagrangian {—R)^ leads only to a field equation of fourth-order for 

D > 2. The lagrangian RU^R with A; > gives a field equation of higher 
than fourth-order. For this case we get 

F = RU^R , G = 2U^R (3.23) 

and no solubility condition for D > 2. For the 00-component of the field 
equation we need 

Fa = U'^-^R (3.24) 

and 

G = 2n''R (3.25) 
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and get 

= d'^R (2R0Q - ^R^ + 2naU^R,Q + 



+ Y.{n'-^R){D^R) - -(rf-^i?),o(n^-^i?),o . (3.26) 

A=l ^ 

The ansatz (3.16) leads to (3.17) and 

n{n^R) = (□'=i?),oo + n(n'=i?),o . (3.27) 
We get the hnearized field equation from (3.26) 

U^R = {U^R)f, . (3.28) 

For A; = 1 we have 

= + + (n^ - n - l)/3 + (-n^ - n)/3 (3.29) 

with the characteristic polynomial 

P(t) = a;^ + 2nr' + (n^ - n - \)x^ + (-n^ - n)a; (3.30) 

possessing the roots = 1, a;2 = 0, = —n and 0:4 = — n — 1. We get the 
solutions 

/5(t) = ci + C2e* + cae-"* + C4e-("+^)* (3.31) 

and 

a{t) = cit + C2e* + cae""* + C4e-("+^)* (3.32) 

and 

lim ^ = 00 (3.33) 
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unless C2 = 0. For the lagrangian RU^R the D-dimensional de Sitter space- 
time is not an attr actor solution of the field equation. The formula 

= (□'=+ii?),o - n'^^^R - (n'i^.o - □'i?),oo + n((n'=i?),o - □'i?),o (3.34) 

for the linearized field equation for A; + 1 leads to the recursive formula for 
the characteristic polynomial: 

characteristic characteristic 

= ■ X ■ {x -\-n) . 

polynomial for k + 1 polynomial for k 
The characteristic polynomial for k has the roots: 





1 


simple 







k-fold 




—n 


k-fold 


Xi — 


-n-l 


simple 



(3.35) 



We get the solutions 

f5{t) = S{t) + T(t)e-"* + cie* + Cac^""-^)* (3.36) 

and 

a{t) = S{t) + r(i)e-"* + cie* + cse^"""^)* (3.37) 

with S, T, T polynomials at most k-th degree and S polynomial at most 
k + 1-th degree. For most of all solutions we get 

lim ^ = oo (3.38) 

and therefore, the de Sitter space-time is not an attractor solution for the 
field equation derived from the lagrangian Rn'^R. 
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These results have shown, that for the lagrangian RU^R with A; > 1 the 
de Sitter space-time is not an attractor solution. The lagrangian (—R)^ 
gives only a fourth-order differential equation. We will try to answer the 
following question: 

Are there generalized lagrangians so, that the de Sitter space-time is an 
attractor solution of the field equation? 
First we make the ansatz 

m 

F^Yl CkRO^R with Cmy^O . (3.39) 

k=l 

In this case the de Sitter space-time is not an attractor solution, because for 
each term there one gets as a root of the characteristic polynomial of the 
linearized field equation. 
Now we make the ansatz 



F = co{-R)^ + CkR^^R with c„ . (3.40) 
for the generalized lagrangian. One gets the characteristic polynomial 



P{x) = x{x + n) 



Co + XI CkX^~\^ + n)''~\x -l){x + n + l) 

k=l 



(3.41) 



for the linearized field equation. The solutions xi = and X2 = —n do not 
depend on the coefficients q of the lagrangian. It is sufficient to look for the 
roots of the polynomial 

m 

P{x) = Co + X CkX^~\x + n)''-\x - l)(a; + n + 1) . (3.42) 

k=l 
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If all solutions of this polynomial have negative real part, then the de Sitter 
space-time is an attractor solution for the field equation. The transformation 

z = x'^ + nx+ — (3.43) 

gives 
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P{x) - Q{z) = 



fe-i 



n \ / n 



k=l+l \ 



2\ fc-'-l 



J-l> 



= do + diz + ... + dmz"" . 

Now let 

da = 1 



(3.44) 



«.* = (-*)""'[(?:;)(-^HY)(f+»+i) 



I = 0, ... ,171 

k — 1, ... ,171 
I < k < m 
else . 



(3.45) 



This gives the equation 
do 



y dm J 



A 



( \ 

Co 



with A regular . 



(3.46) 



The roots of P{x) have a negative real part iff the roots of Q{z) are from the 
set 

I T)^ / T7 2 I 

(3.47) 



Tl^ / Hp' 

M := + iy : X > — /\\y\ < n\l x — — 



If the roots Zk of the polynomial Q{z) are elements of M, then the coefficients 
dk are determined by 



Q{z) = Y,dkZ^ = IK^-^fe)- 



(3.48) 



fc=0 



k=l 
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The coefficients 



Co 



y c^ J 



do 



(3.49) 



belong to a lagrangian, that gives a field equation with a de Sitter at- 
tr actor solution. The above considerations have shown that for every m 
there exists an example for coefficients Cfe, so that the de Sitter space-time 
is an attractor solution for the field equation derived from the lagrangian 

co(-i?)^ + E CkRU^R with Cm ^ 0. 

fc=l 

It turned out that all the variants of the definition of an asymptotic de 
Sitter solution given in subsection 2.1. lead to the same class of solutions. 

For the 6th-order case we can summarize as follows: 
Theorem 1: Let L ^ Ej^ + ciRUR and Le ^ R - ^-^ L with length 
/ > 0. Then the following statements are equivalent. 

1. The Newtonian limit of is well-behaved, and the potential 4> consists 
of terms ^e""'" with a > only. 

2. The de Sitter space-time with i7 = y is an attractor solution for L in the 
set of spatially flat Friedman models, and this can already be seen from the 
linearized field equation. 

3. Ci > and the graceful exit problem is solved for the quasi de Sitter phase 

H < l/l of Le. 

4. l'^ — 1'^ + ll, Pci — has a solution with < Zq < h. 

5. < ci < f . 
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For the proof of 3. one needs eq. (18) of the first of papers [64] which reads 
in our notation H{1 — AciH"^) — —1/6P showing that ii" < at the quasi de 
Sitter stage. 

Theorem 2: Let L and Le as in theorem 1. Then are equivalent: 

1. The Newtonian hmit of is well-behaved, for the potential 4> we allow 
i and terms like ^^e'""^ with a > and a polynomial P. 

2. The de Sitter space-time with H = j cannot be ruled out to be an 
attractor solution for L in the set of spatially flat Friedman models if one 
considers the linearized field equation only. 3. Le is tachyonic-free. 

4. P + ll, fci = lUl has a solution with < < ^i- 5. < Ci < ^. 

Of course, it would be interesting what happens in the region where the 
linearized equation does not suffice to decide; one should even not try to 
answer this question without a computer algebra system. 

4 Higher-order gravity and diagonalizability 
of Bianchi models 

A Bianchi model can always be written as 

ds'' = de - gap{t)a"al' (4.1) 

where Qafj is positive definite and cr" are the characterizing one-forms. It 
holds 

da^^-\clpa'^^a^ (4.2) 
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with structure constants C^^ of the corresponding Bianchi type. It belongs 
to class A if Cf^ = 0. The abelian group (Bianchi type I) and the rotation 
group (Bianchi type IX) both belong to class A. 

In most cases, the Qaf} are written in diagonal form; it is a non-trivial 
problem to decide under which circumstances this can be done without loss 
of generality. 

For Einstein's theory, this problem is solved in [46]. One of its results 
read: 

If a Bianchi model of class A (except types I and II) has a diagonal 
energy-momentum tensor, then the metric gapit) can be chosen in diagonal 
form. Here, the energy-momentum tensor is called diagonal, if it is diagonal 
in the basis {dt, a^, a^, a^). 

This result rests of course on Einstein's theory and cannot be directly 
applied to higher-order gravity. 

For fourth-order gravity following from a Lagrangian L — f{R) considered 
in an interval of i?- values where 

df <Pf , 
— ■ — - 4 
dR dR^ ^ 

one can do the following: The application of the conformal equivalence theo- 
rem is possible, the conformal factor depends on t only, so the diagonal form 
of metric (3.1) does not change. The conformal picture gives Einstein's the- 
ory with a minimally coupled scalar field as source; the energy-momentum 
tensor is automatically diagonal. So, in this class of fourth-order theories of 
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gravity, we can apply the above cited theorem of MacCallum et al. 

As example we formulate: All solutions of Bianchi type IX of fourth-order 
gravity following from L — B? considered in a region where R ^ can be 
written in diagonal form. 

Consequently, the ansatz used in [27] by Barrow and Sirousse-Zia for this 
problem is already the most general one, cf. [30] Spindel. 

For fourth-order gravity of a more complicated structure, however, things 
are more involved; example: Let 

with ab 0. Then there exist Bianchi type IX models which cannot be 
written in diagonal form. (This is a non-trivial statement.) 

To understand the difference between the cases 6 = and 6 7^ it proves 
useful to perform the analysis independently of the above cited papers [46] . 
For simplicity, we restrict to Bianchi type I. Then the internal metric of the 
hypersurface [t = 0] is fiat and we can choose as initial value gapi^) = 
Spatial rotations do not change this equation, and we can take advantage of 
them to diagonalize the second fundamental form ^ga/3{0). 

First case: b — 0. As additional initial conditions one has only R{0) and 
^i?(0). The field equation ensures gapit) to remain diagonal for all times. 

Second case: 6 7^ 0. Then one has further initial data -^gapiQ)- In the 
generic case, they cannot be brought to diagonal form simultaneously with 
^5*0/3 (0). This excludes a diagonal form of the whole solution. (To complete 
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the proof, one has of course to check that these initial data are not in con- 
tradiction to the constraint equations.) This case has the following relation 
to the above cited theorem [46]: Just for this case b ^ 0, the conformal re- 
lation to Einstein's theory breaks down, and if one tries to re-interpret the 
variational derivative of Cij^iC'^^^'' as energy-momentum tensor then it turns 
out to be non-diagonal generically, and the theorem cannot be applied. 

For higher-order gravity, the situation becomes even more involved. For 
a special class of theories, however, the diagonalizability condition is exactly 
the same as in Einstein's theory: If L = R + Y.k=o dkRO^R, {am 7^ 0) then in 
a region where 2L ^ R the Cauchy data are the data of General Relativity, 
i?(0), and the first 2m -\- 1 temporal derivatives of i? at i = 0. All terms with 
the higher derivatives behave as an energy-momentum tensor in diagonal 
form, and so the classical theorem applies. [Let us comment on the restriction 
2L ^ R supposed above: Eqs. (3.5, 3.6) show that Fq — G — Q represents a 
singular point of the differential equation (3.3); and for the lagrangian given 
here G = ^ — 1. For fourth-order gravity defined by a non-linear lagrangian 
L{R) one has G — ^ and G — defines the critical value of the curvature 
scalar.] 
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5 Structural stability of fourth-order cosmo- 
logical models 

In [23] , Coley and Tavakol discuss cosmological models from the point of view 
of structural stability; the notion for the contrary of it is fragility Structural 
stability is a more general but less strictly defined notion than the usual 
stability So, its concrete meaning has context-dependently to be specified. 

1. Example: The Einstein universe (a closed Friedman model of constant 
world radius in General Relativity with positive cosmological term ~ A and 
incoherent matter as source) is unstable with respect to the initial data: A 
non- vanishing but arbitrarily small initial Hubble parameter gives rise to a 
singularity This property ruled out the Einstein universe as describing our 
real world. It should be emphasized that this is in coincidence with the 
observational result that our universe is not static, but that this theoreti- 
cal stability analysis ruled out the Einstein universe independently of the 
observational result. 

Structural stability represents stability not only with respect to a small 
perturbation in the initial data, but a small change in the corresponding type 
of matter and field equations. In most of the specifications one requires that 
by a small change of conditions the qualitative (or topological) properties 
of the system remain unchanged. Concerning field equations, Coley and 
Tavakol [23] concentrate on Lagrangians L — f{R) for the gravitational field: 
For linear functions / one gets General Relativity, for non-linear ones fourth- 
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order gravity. Because of the change in the order of the differential equation 
the question concerning the robustness of General Relativity is a non-trivial 
one. Before we follow this line we present some more or less trivial examples 
from General Relativity to clearify what is meant. 

2. Example: The spatially flat Friedman model with incoherent matter 
(dust) but A = (Einstein-de Sitter model) has a scale factor a ~ t^^^ for 
synchronized time t. A small change of the initial data only changes the 
proportionality factor, so this is stable. However, if we consider this model 
within the class of all Friedman models, then it represents just the bifurcating 
point between the ever-expanding open and the recoUapsing closed models. 
In this sense, the Einstein-de Sitter model is a fragile one. 

3. Example: Again we consider the Einstein-de Sitter model within the 
class of all spatially fiat Friedman models. We impose new structure by 
allowing a new contribution to the energy-momentum tensor in form of ra- 
diation not interacting with the dust. During expansion, the energy density 
of the radiation falls ~ a~^, and of the dust only ~ . The radiation be- 
comes asymptotically negligible, and, asymptotically for large values one 
gets approximately a ~ t^/^. In this sense, the Einstein-de Sitter model is 
structurally stable. 

4. Example: Now we invert the point of view from the third example. 
We start from a spatially flat Friedman filled with radiation. Then one has 
a ~ i^/^. We impose new structure by adding an arbitrarily small amount of 
non-interacting dust. As in the previous example, we get asymptotically a ~ 

36 



t^l^. In this sense, the spatially flat Friedman radiation model is structurally 
unstable. 

Let us now come to the consideration of Coley and Tavakol concerning 
structural stability of fourth-order gravity models. They consider pertur- 
bations of Friedman's radiation model within fourth-order gravity. For the 
non-tachyonic case, they get as result that the R^-ierva. gives rise to an in- 
stability. It is known for a long time, that asymptotically 

the i?^-term gives rise to damped oscillations which behave as dust in the 
mean. So, the structural instability considered there is exactly the same as 
in the 4. example above and not a special feature of the fourth-order term. 

Analogously they consider the quasi-de Sitter stage (Starobinsky infla- 
tion) and get its stability for the non-tachyonic case L = R + aE?. 

Remark: Coley, Tavakol [23] use the notion "topological almost all" in the 
sense of " countable intersection of open dense subsets" . One must be careful 
in applying this notion, especially, if one is tempted to mix it with the notion 
"almost all" in measure theory. A remarkable example shall underline this 
warning: Let / = [0, 1] be the closed interval with the usual probability 
measure ii. Let {rn\n e A^} C / be a countable dense subset of /. For each 
natural m, let 

neN 

and A — flmeAr ^m- (Here, ]x, y[ denotes the open interval.) Each is open 
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and dense in I. For all values m, 

Hence, = 0. So, A contains topologically almost all points of /, but 

there is zero probability to meet an element of it. 

Next, Coley and Tavakol consider structural stability of Starobinsky in- 
flation L = R + aE? with respect to addition of the cubic term hB?. For 
L — alone one gets polar inflation a ~ t~^^] considered in the region 
t < this is expanding with hh~'^ = ^. The term with b does not alter the 
order of the differential equation, and so one expects a continuous change of 
the properties. In fact, for small values |6| one has Starobinsky inflation as 
transient attr actor, with increasing |6| one gets a smaller basin of attraction, 
and for |6| ^ one needs fine-tuned initial conditions. 

A more drastic change of structure is to be expected if we consider struc- 
tural stabihty with respect to the addition of terms like RnR. 

6 Discussion 

Sudarsky [66] proves the no-hair theorem (in the version that there are no 
non-trivial black holes with regular horizon) for the Einstein-Higgs theory. 
We have deduced a cosmic no hair theorem on a quite different footing as 
follows (the more detailed formulation is given at the end of set. 3) 
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Theorem: Let L = B? + RUR and 

Le = ~ f' "^ith length / > and arbitrary real e. Then the 

following statements are equivalent. 

1. The Newtonian hmit of is well-behaved. 

2. The de Sitter space-time with = j is an attr actor solution for L. 

3. e > and the graceful exit problem is solved for the quasi de Sitter phase 
H < l/l of Le. 

4. P — 1^ + ll, el^ — l^ll has a solution with < lo < h- 

5. < e < i. 

Prom the first glance this theorem is contrary to the results of refs. [64 
- 66]. But one should remember that in refs. [64 - 66] the question had 
been considered whether the sixth-order terms can typically lead to double 
inflation. The answer was: Double inflation (one period from the i?^-term, 
the other one from the ROR-term) requires a flne-tuning of initial conditions. 
Here we have shown: The results of the Starobinsky model (e = in the 
present notation) are structurally stable with respect to the addition of a 
sixth-order term ~ eROR^ where < e < |. The duration of the transient 
quasi de Sitter phase becomes reduced by a factor ~ (1 — 4e) only. 

Further we have shown: For 

L ^ R^ + CiRnR + C2ROOR, C2 7^0 

and the usual case n — 3 the de Sitter space-time with = 1 is an attractor 
solution in the set of spatially flat Friedman models if and only if the following 
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inequalities are fulfilled: 

1 1 

< Ci < -, < C2 < — 

4 16 

and 

ci > -13c2 + ^4c2 + 225ci 

This represents an open region in the Ci — C2-plane whose boundary contains 
the origin; and for the other boundary points the linearized equation does 
not suffice to decide the attractor property. This situation shall be called 
" semi-attractor" for simplicity. In contrary to the 6th-order case, here we 
do not have a one-to-one correspondence, but a non-void open intersection 
with that parameter set having the Newtonian limit for Le well-behaved. 

To find out, whether another de Sitter space-time with an arbitrary Hub- 
ble parameter i7 > is an attractor solution for the eighth-order field equa- 
tion following from the above Lagrangian, one should remember that H has 
the physical dimension of an inverted time, Ci a time squared, C2 a time to 
power 4. So, we have to replace ci by ciH^ and C2 by C2i?^ in the above dimen- 
sionless inequalities to get the correct conditions. Example: < ciH'^ < |. 
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